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The beneficial effects claimed for whirl control devices 
demonstrate that the dynamic behaviour of rotors is Influenced by 
the fluid whirl in shaft and balance drum seals. The present paper 
reports results from two series of experiments, the first on the 
factors affecting the whirl at the seal inlet, and the second on 
the variation of whirl velocity along the seal. In both cases the 
LDA measurement technique required the clearance between the fixed 
and rotating parts of the models to be substantially greater than 
occurs in real machines, but the results are indicative 
nevertheless. Experimental and theoretical results are given for 
the radial distribution of whirl velocity in the gap between 
impeller shroud and pump casing. Results of tests with modified 
stator surfaces are also shown. This work leads naturally Into the 
second series of experiments where some preliminary measurements of 
velocity distribution in the clearance between a fixed stator and a 
rotating shaft are reported for a range of inlet whirl conditions. 



NOTATION 

Co Flow parameter Q/flR 2 h 

f Friction factor 

h Clearance gap size 

k Constant 

Q Leakage flow towards axis 

R Impeller radius 

r Radial coordinate 

t Time 

U Mean axial velocity in annulus 

v Tangential velocity 

vf Friction velocity 

z Axial coordinate 

e Eddy kinematic viscosity 

v Kinematic viscosity 

p Fluid density 

t Shear stress 

Q Angular velocity of impeller 

« Angular velocity of fluid 


Suffices 


1 Impeller shroud surface 

2 Surface opposite to impeller shroud 

d Edge of balance drum 

R Edge of impeller 

s Shaft surface 


* The work reported here was partially supported by Weir Pumps Ltd. 
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INTRODUCTION 


The part played by fluid whirl In determining the stiffness 
characteristics of shaft and balance drum seals Is widely 
acknowledged (1,2, 3, 4). The present report applies particularly to 
the situation at the high pressure end of a multi-stage centrifugal 
pump, where a small fraction of the flow leaving the last Impeller 
leaks away through the shaft or balance drum seal. The diameters of 
shaft, drum and Impeller are typically In the ratios 1:2:3 and the 
leakage path is shown schematically In Fig. 1. 

For both experimental and analytical purposes the arrangement 
shown In Fig. 1 can be conveniently Investigated in two obvious 
sections, firstly the whirling radial Inflow between the impeller 
shroud and the casing, and secondly the whirling axial outflow 
through the seal. The stiffness characteristics are determined by 
events in the seal, but the radial flow in the gap between shroud 
and casing is important because it influences the whirl of the 
leakage flow at entry to the annular seal. A number of whirl 
control devices have been found to have beneficial effects on rotor 
stability (5,6,7), thus demonstrating the inaccuracy of earlier 
theories in which the flow entering the seal was supposed to 
acquire instantly a whirl velocity independent of the entry value. 

When a rotor surrounded by fluid in a small clearance gap moves 
off-centre, it experiences a force which depends on both the 
position and the translation velocity of the shaft centre, and 
which has both radial and tangential components. In particular. 
Black et al (1) showed that a radial displacement produced a 
tangential force component dependent on the fluid whirl velocity. 
The interaction between rotor and fluid is clearly extremely 
complex, and one way to simplify the problem of computing rotor 
motion is to make some plausible assumption about the variation of 
whirl along the seal. The force on the shaft is then obtained by 
integration of the contributions from all the axial elements of the 
seal, each depending on the assumed local whirl velocity of the 
fluid. The proposal of Black and co-workers for an asymptotic 
approach towards a limiting value of whirl appears to be a useful 
development beyond the discredited step change at entry mentioned 
above . 

Obvious difficulties arise with this approach when the rotor 
and seal do not have a concentric equilibrium position and fluid 
particles experience different conditions on their passage through 
the seal. In some cases the short bearing assumption may be invoked 
to Justify treating the whirl velocity as independent of angular 
position, but experimental evidence is hard to find. 

IMPELLER/CASING FLOW 

The angular velocity of the flow leaving the impeller and entering 
the gap between shroud and casing is typically about half that of 
the shaft, although it can be much higher under low flow 
conditions. The clearance gap is several millimetres wide and the 
situation may be regarded as a steady axl-symmetric flow provided 
that shaft misalignment and bending do not cause significant 
angular or temporal variations in the clearance. The flow will be 
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The whirl profile occurring in any real case is modified by 
further effects which theoretical work may need to account for _ Two 
mechanisms tend to reduce variations of angular velocit^The first 
operates by the action of viscosity which produces tangential sh 

stress between adjacent annular elements of fluid Lg u ®ualW 
is the turbulent momentum exchange between elements usually 
described as Reynolds shear stress. Viscous stress is neglected 1 

the P r ®® C ^ t h ^ n ^ f S e 1 ( ft arises from the secondary flows adjacent to 
the boundaries. The whirl approaches 0 in the thin fluid layer next 
¥n fho 1 mnpl ler shroud and, as this is normally faster than the 
main flow at the same radius, the faster moving fluic d to f ^* 

outwards. Similarly the slow-moving fluid next to _ the fl d 
boundary tends to move Inwards. In cases where the leakag 
zero or small, these secondary flows can be an important means of 
transporting angular momentum in the radial directi . 
transport is always outwards, whatever the gradient of e ma ^ 
SI Profile The effect of secondary flows will be small at 
higher leakage rates except where the net radial flow Is zero, in 
the region : 


r < r < r . 
s d 


RADIAL FLOW RIG 

Fig. 3 shows the arrangement of the test rig used to investigate 
j • _ i inwards flow The vertical shaft carries two discs, one 
200mm diameter representing the Impeller shroud and a 
r^nresent ine the end of the balance drum. The inlet to the seal is 
represented by the clearance between the edge of the small disc and 
an P annular plate representing the pump casing. Spacers onthesh ^ 
can be changed to vary the effective shaft diameter and the gap 
between the 8 dlscs. The test liquid is water and the valves on two 
Dumped loops are used to control the simulated leakage flow and its 
uMrl asifenters the gap at the edge of the large disc. The whirl 
is induced by the outflow from a whirl chamber with a number of 

tangentiali inlets. thrQUgh the sides of the glass-walled 

containing tank. Velocity measurements are made by forward scatter 
laser doppler anemometry with the optical axis of the L y 
passing horizontally through the gap between the discs. The beam 
intersection can be moved on three axes, but direct measurement of 
[he whirl component of velocity by placing the optical axis along a 
disc diameter is unfortunately obstructed by the shaft. Velocity 
components have to be computed from measurements at off diameter 
positions of the optical axis. 

RADIAL FLOW RESULTS 

Typical velocity distributions of whirl across the gap width are 
shown in Fig. 4, where the angular velocity is virtual y un 
except for thin layers on the walls. These findings provide 
justification for restricting most further measurements to the 
mid-plane of the gap and also for a one-dlmenslonal analysis. The 
Reynolds numbers given on the figures are values of OR /«■ 

Whirl profiles for a series of tests in which only th 
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rate Q was varied are shown in Fig. 5. Increasing the flow rate 
alters the balance between frictional and inertial effects by 
increasing the rate at which angular momentum is convected into the 
gap. The result is a general increase in whirl velocity. 

Fig- 6 shows the effect on the whirl profile of modifying the 
stator surface with other controlled variables held constant. The 
highest line is the standard nominally smooth surface and the 
dashed line shows the effect of a replacement surface with blind 
holes 3mm dia x 3mm deep on a 7mm pitch. The other two curves were 
obtained for stators with twenty-four radial slots 6mm wide 
extending out to 0.76 x impeller radius. It may be seen that these 
last results diverge significantly from the smooth surface result 
only downstream of this radius. 


RADIAL FLOW THEORY 

A one-dimensional computer model was written to solve the whirl 
profile, based on the angular momentum equation for an annular 
element of axial height h and infinitesimal radial extent: 

r 3 -aw/at = Q/2irh-a(r 2 w)/ar + r 2 (x 2 -xi)/ph - a(r 2 T/p)/3r. (6) 

Shear stresses ti and X 2 are applied to the -flow by the opposing 
walls: 

Ti/p = fi -r 2 (ni-u) (ni~«)/2 (7) 

and X2/p = f2T 2 (w-n2)(w~02)/2 (8) 

where ~ denotes the positive difference between two quantities. 
The Reynolds stress x, acting on the interface between adjacent 
elements can be written as the product of an effective viscosity 
and the velocity gradient: 

x = -pe r-dw/ar. ( 9 ) 

In reality c will vary across the gap in some way such as: 

e = k(vn + vf 2 ) - z(h— z)/h (10) 

which at least ensures that it is zero on both walls and gives a 
mean value for use in eq. (9) as: 

e = (kh)(vn + vf2)/6. ( 11 ) 

Each friction velocity can be related to the friction factor on the 
corresponding wall and the relative velocity. Thus, considering 
rotational velocities only: 

vn 2 = ( f i/2 ) [ (fii~w)r] 2 and vf 2 2 = (f2/2) [(w~fi 2 )r J 2 (12) 
which allows c to be written as: 

e = kh/6*fir (f/2) 0 ' 5 where f 0-5 n = f i 0-5 (ni~u) + f 2 0,5 (w~ 02 ). (13) 

The friction factors may be treated as constants or as functions of 
a local Reynolds number. 

Using equation (6) to follow the progress of the transient from 
an arbitrary initial whirl profile can be made the basis of a 
numerical solution method. The accuracy of following the transient 
does not affect the steady-state solution. The equation can be 
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conveniently made dimensionless by using R and £1 as reference 
quantities so that a>/tl becomes o>, r/R becomes r, and lit becomes the 
shaft rotation angle from the start of the transient. If 
differentiation Is denoted by a dash for radius and a dot for 
angle, equation (6) can be written In dimensionless form as: 

r 3 u = (Q/ZnflR^) ■ (r 2 w)' - r 2 * (t2-ti J/ptl^h ~(r 2 • T/pfl^ 2 ) ' . (14) 

The upper limit on the angular step required to avoid numerical 
Instability in the discretlsed form of this equation with 20-30 
radial steps is typically 5-30 degrees. There Is a rapid approach 
to within a few per cent of the steady state solution in the region 
where r>rD, but further progress is then limited by slow 
convergence in the diffusion-controlled region near the shaft where 
r<ro. 

Figs. (7,8) are whirl predictions typical of actual machines and 
the test rig, which show similar behaviour to the experimental 
results. A feature observed in some of the test flows is a region 
of positive w' near the edge of the balance drum, even though the 
fluid whirl Is below shaft speed. This cannot be accounted for In 
the one-dimensional model and must be attributed to the 
two-dimensional character of the flow as it turns Into the annular 
seal. 


ANNULUS FLOW 

As previously discussed, some rotor dynamicists have used an 
assumed whirl profile along the seal clearance, and one the main 
objectives of the present work was to Investigate the form of the 
actual profile. Initial experiments were therefore carried out In a 
concentric system with and without shaft rotation, and with and 
without Inlet whirl. The assumption of a whirl profile with no 
angular variation can be strictly true only in the concentric case, 
but It has been extended to eccentric cases where the short bearing 
theory could normally be justified. The second aim of the work was 
to test the validity of this extension by measuring the tangential 
and axial velocity components at a large number of points In the 
eccentric annulus. 


ANNULUS TEST RIG 

Fig. 9 shows the test rig used to study the behaviour of whirling 
flow in an annulus. The arrangements for controlling the leakage 
flow rate and the inlet whirl are broadly similar to those used In 
the radial flow rig, but access to the flow for LDA measurements is 
more troublesome. A narrow window was fitted along the length of 
the plain stator and a precision glass tube mounted in roller 
bearings was used to represent the shaft. It was driven at one end, 
leaving the other end open for an angled mirror to be inserted. 
Thus the optical axis is folded, with the transmission radial and 
the collection axial. A frequency shifter was used to resolve low 
velocities. The rotor diameter is 92mm and measurements were made 
over 140mm axial length. The mean radial clearance is 10mm which is 
clearly very large in comparison with real seals but the Reynolds 
number and the inlet whirl ratio are however comparable with values 
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in actual machines. 

Two series of tests have been carried out, one series with tho 
rotor and stator concentric, and the other 

cases™ S "" re ” ade at SlX angular Positions In the eccentric 




Results from some concentric tests are shown in Figs 10 and 11 
where measures of tangential velocity are plotted against axil 
position along the annulus for fixed flow rates. The data for 

velocity” useful? ^ Shaft statlonar y- 50 the tangential 

[Tflr /h k f ! expressed in terms of the mean axial velocity 

the ratio Tf The 0 fi lnlet W ? irl conditions which are quantified by 
totaT f n tK %-f 1 passing through the swirl inducer Qs to the 

w Qt. The axial Reynolds number (2Uh/i>) is tk 

a"o VeS Th f o 0r ap h pa 8 r:n r t W ?i lrl T ‘ ^°-a%ds “opaS 

certainly X to l °V Curves ls •*“» 

tests ibUUOn WhlCh haS ^ measured ln dTuiT™ in ln S ome e of^the 

294^") ^ 8i " * e ‘ ^h 'was Girled whV’lV ^ V“\l 5 

dS^tL^=‘ lal — ' S» K 


( V W eo )/(u o ” "» )= ej «P(-ax) (15) 

where ^ is the limiting angular velocity as x tends to infinity. 

Taking x as the test section length L, equation (15) can be written 


V <V (1-C)w «o where C = exp(-aL). (i 6 ) 
When data from Fig. 11 is plotted as u/£l against « /Q, the points 

lie closely about a straight line from which estimates of C and u 
can be derived: 00 


w m = o.4 

oo 


ana 


C = 0.42. 


between that the leaka « e flo » 
velocity at the seal ini .f P ™ P CaSing Can result in « whirl 
one-dimensional theory c^n mode? £“ ter “T half Shaft S P^. A 
with empirically determine? friction factors behaVi ° Ur of the flow 
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Different stator surface modifications were tested and radial 
grooves were found to be an effective whirl-control device. 
However, one effect of reducing the whirl of the radial Inflow will 
be to increase the pressure at the seal entry. The leakage will 
therefore tend to increase and the whirl will consequently 
increase. This self-compensating effect of attempts to reduce the 
whirl can be minimised by concentrating efforts to modify the flow 
in the region Just before the flow enters the seal. 

In the concentric annulus flows so far tested, the influence of 
the inlet whirl persisted along the length of the test section and 
could be described quite well by an exponential approach to a 
limiting value. 
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Fig. 2 Types of whirl profile. 
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Fig. 3 Radial flow rig. 
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Fig. 4 Axial variation of whirl 
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Fig. 5 Radial variation of whirl. 
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Fig. 6 Whirl profiles for different stator surfaces. 



Radius (r/R) 

Fig. 7 Theoretical whirl profiles for test rig. 



Fig. 8 Theoretical whirl profiles for pump. 
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Fig. 10 Tangential velocity variation along a fixed annulus. 



Fig. 11 Whirl variation along annulus with rotating shaft. 
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